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Anisotropic pairing interactions mediated by phonons are examined in layer systems. It is shown 
that the screening effects become weaker when the layer spacing increases. Then the anisotropic 
components of the pairing interactions increase with the screening length since the momentum 
dependence changes. As a result, various types of anisotropic superconductivity occur depending on 
the parameter region. For example, p-wave superconductivity occurs when the short-range part of 
Coulomb repulsion is strong and the layer spacing is large. Two kinds of inter-layer pairing may occur 
when the layer spacing is not too large. Although the phonon contribution to the d-wave pairing 
interaction is weaker than the p-wave interaction, it increases with the layer spacing. Relevance of 
the present results to organic superconductors, high-Tc cuprates, and Sr2Ru04 is discussed. 



I. INTRODUCTION 

Anisotropy of superconducting order parameter and 
the mechanism of pairing interactions in layered super- 
conductors are recent subjects of much interest. In par- 
ticular, high-Tc cuprate superconductors, organic super- 
conductors, and Sr2Ru04 compound have been studied 
by many authors. 

There are some evidences that the order parameter 
has line nodes on Fermi-surface in high-Tc cuprates. For 
example, an experiment and a theory on Josephson junc- 
tion gave an evidence of "d-wave" order parameter in 
a cuprate superconductor [0j2j. Linear temperature de- 
pendence of the penetration depth was observed at low 
temperatures |3|. 

On the other hand, superconductivity in Sr2Ru04 is 
considered to be due to spin triplet pairing according 
to the results of a Knight shift measurement Q and 
a ^iSR experiment In the organic superconductor 
(TMTSF)2PF6, a Knight shift measurement seems to 
support spin triplet pairing 1^,0. 

On mechanism of the high transition temperature of 
the cuprates, pairing interactions of magnetic origin, 
such as exchange of spin fluctuations and a super- 
exchange interaction between nearest neighbor spins, 
have been discussed by many authors, because of proxim- 
ity to the antiferromagnetic phase. However, experimen- 
tal results of isotope effect suggest that there are contri- 
butions to the superconductivity from phonon-mediated 
interactions in many high-Tc cuprates |Tc| [T^ . Absolute 
values of shifts of Tc are very large (0.2K ~ 0.7K), but 
isotope effect exponents a are small because of the high 
transition temperature. 

Abrikosov proposed a theory based on weak screen- 
ing of Coulomb interactions and phonon-mediated pair- 
ing interactions in which anisotropic s-wave order pa- 
rameter was obtained [0. In the presence of on-site 
Coulomb repulsion, extremely anisotropic s-wave order 
parameter with nodes was obtained |lq]. Bouvier and 



Bok also calculated an order parameter explicitly, and 
obtained anisotropic s-wave in the same model 1 1^ . Re- 
cently, it has been shown that d-wave superconductivity 
is reproduced in a similar model with antiferromagnetic 
fluctuations 

We proposed in our previous paper that triplet 
pairing superconductivity can be induced by phonon- 
mediated interactions in ferromagnetic compounds, 
where singlet pairing is suppressed by Pauli paramag- 
netic effect. 

The origin of the anisotropic components of pairing 
interactions mediated by phonons is briefly explained as 
follows. The screening effect limits electron-ion interac- 
tions within a range of the order of the screening length. 
Since the pairing interactions mediated by phonons are 
obtained by a second order perturbation of the electron- 
ion interactions, they also have a range of the same order. 
For example, the screening effect is taken into account 
as vertex corrections within diagramatic technique [ p3| . 
When the screening length increases, the interactions 
are more localized in the momentum space. Hence, the 
anisotropic components of the interactions increase with 
the screening length. 

In this paper, we examine layered superconductors 
with the phonon-mediated pairing interactions, extend- 
ing it to systems with large layer spacing. The layered 
structure modifles the screening length and the pairing 
interactions signiflcantly. It is shown that anisotropic 
components of the pairing interactions increase with the 
layer spacing. We argue from this result that some as- 
pects of the layered superconductors, such as Sr2Ru04, 
organics, and cuprates, can be explained in the present 
model. 

We also study an effect of anisotropy of density of 
states in square lattice systems. Although the effect of 
the anisotropy must be most remarkable when the Fermi- 
surface is near the van Hove singularities, we consider 
a system not necessarily near the van Hove singularity 
but a system with the density of states anisotropy within 
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the layers. Enhancement of the transition temperature 
due to the van Hove singularities has been discussed by 
many authors. In particular, it has often been discussed 
that the effect is efficient especially for the d^2_y2 pair- 
ing p^ , p5| ,p|, p6y20|j2l[ | . In this paper, we concentrate on 
the anisotropy but do not discuss the value of the density 
of states. It is shown that superconductivity is enhanced 
by the anisotropy for pairing of d^2_y2 symmetry, but 
not for those of dxy, Px, and Py symmetries. 

In section II, we define the model of the pairing in- 
teractions mediated by phonons. We derive expressions 
of the coupling constants for various types of anisotropic 
superconductivity. In section III, we examine the depen- 
dence on the layer spacing of the screening length and 
the pairing interactions. In section IV, we consider a 
situation in which inter-layer coupling is of the order of 
intra-layer coupling. In section V, we examine an effect 
of the anisotropy in the electron dispersion in square lat- 
tice systems. Section VI is devoted to discussion and 
summary. 



eq. (2) in the form 



y(q) = - 
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qiiP + 9s2 



cos q^b (3) 



for layer systems, where qy is the momentum element in 
the layers. Here we have truncated the interaction at the 
nearest layers. The first and second terms correspond to 
the intra-layer and inter-layer interactions, respectively. 
The parameter q'^ is the inverse of the range of the inter- 
actions between electrons on the nearest layers in the x 
and y-directions. 

The gap equation of supercoductivity is written as 



where 



Wik') - 2T 



2E{k') 



(4) 
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II. SCREENING EFFECT AND PAIRING 
INTERACTIONS 

First, we introduce a model of pairing interactions. 
Abrikosov examined an effective pairing interaction me- 
diated by phonons of the form 



^(q) = 5(- 



J' 



(1) 



with q = Iql and — l^^, where Ig denotes the screening 
length |17| , |l8|| . A similar form corresponding to rt = 1 
is obtained by taking into account the screening effect 
in electron-phonon interactions as explained in a text 
book If we put n = 1 for simpHcity and ^k— Ck+q = 
for the electrons near the Fermi-surface in eq. (1) accord- 
ing to Abrikosov , we obtain a simplified form 



V'(q) 



gqs 



We examined this model in three dimensions in our pre- 
vious paper p2| . 

In this paper, we consider the layer systems, and define 
lattice constants a within the layers and b between the 
layers. We take x and y-axes in the direction of the lat- 
tice vectors within the layers, and z-axis perpendicular 
to the layers. In eq. (2), the range of the interaction is 
Is/b layers in the z-direction, while it is Is /a sites in the 
layers. Thus, the anisotropy due to the layered structure 
is partially taken into account for the difference of the 
lattice constants b and a. This model is qualitatively ap- 
propriate for long wave length such as A ~ q^^ 3> a, b, Ig. 
However, when 5 3> a, the discrete layered structure in 
the inter-layer direction must be taken into account for 



shorter wave length \ ^ q 



b. Therefore we extend 



with £'(k) = y/el_ + [A(k)]2 and N the number of lattice 
sites. 

We put the gap function 



A(k) = A||(k||)r;(fc, 



(6) 



where ky = (k^, ky) and ri{kz) is a normalized function of 
the momentum component fc^. From eq. (3), the solution 
of the gap equation (4) at T = Tc has a form with ri(kz) = 
1, v^cosfczfo, or \/2sinfcz6. Then eq. (4) is written as 



k' 



where N\\ denotes the number of sites in a layer, and 
V^(k||,k'|) denotes the averaged pairing interaction de- 
fined by 



(2) VO^h^\) 



r/b 



TT/b 

dfcz 

— 7r/b J—7z/b 



dfc;r;(fc,)T/(k,k')r;(fc;). 



(8) 

Here we assume that the dispersion in the z-direction can 
be neglected in ek in the gap equation. 

We consider cylindrically symmetic Fermi-surface from 
now on. Hence we put |k||| — |k||| — k-p in the pairing 
interactions eq. (8) and obtain 

v{^~^') ^ v{kM = — ^„^"r^^ (9) 



a — cos((y9 — (p') ' 



with 



a = 1 



2kp 



(10) 



for ri{kz) = 1. On the other hand, for the order parame- 
ters with ri{kz) — ^/2cos kzb and ri{kz) = v^sinfc^fo, the 
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expression for V{(p — (p') is obtained by replacing g and 
a with g' and a' = 1 + q'^ /2fc|, respectively, in eq. (9). 
We expand the averaged interaction V{ip — ip') as 

C30 

V{(fi-ip') = ^ Vmnmlmi^ - 

oo 

m=0 

(11) 

and the gap function A||((/?) = A||(kj|) as 

oo 

m (12) 

m=0 

where 

I -7^ fo?) = cosfmw") 

(13) 



lm{^) = n„cos(m<^) 
7m (<p) = nmSin(m</9), 



with normalization factors 

1 for m = 
\/2 for m 7^ 0. 

The expansion factor Vm is calculated by 



Vm = — r^lm{0)V{e). 
nm Jo ^'^ 



(14) 



(15) 



It is easy to perform the integration in eq. (15). For 
V{kz) = 1, we obtain dimensionless coupling constants 



An 



(16) 



Then the superconducting transition temperature Tc is 
obtained by 



Tc = 1.13a;De 



-l/Xn 



(17) 



with Am = -VmN(O) from cq. (4), where iV(0) is the 
density of states per site of a given spin. 

For r]{kz) = \/2 cos kzb and \/2 sin k^b we obtain a sim- 
ilar dimensionless coupling constant as 



= ^5'iV(0) [a'-V^iy 



(18) 



for nearest neighbor layer pairings. The expression for 
Tc is the same as eq. (17). 

Here, we note that a contribution from the short-range 
part of the Coulomb repulsion must be subtracted from 
Ao obtained above. For example in the tight binding 
model, the on-site Coulomb energy is estimated by 



U = JJ d^rd^r'\w{r)\^' 



47reo|r — r 



Hk(r')r, (19) 



where w(r) is the Wannicr function. It is obvious that 
the energy U is not included in our interaction energy 
eq. (3), since eq. (19) depends on the profile of the Wan- 
nier function. Equation (3) describes the behaviors of 
pairing interactions of longer wave length, while the en- 
ergy U in eq. (19) is characterized by the local states of 
electrons on each lattice site. 

Therefore we must consider the on-site Coulomb re- 
pulsion in addition to the pairing interaction of cq. (3). 
However, it reduces only the intra-layer s-wave pairing 
interaction but not the other anisotropic pairing inter- 
actions because of the symmetry. We define a param- 
eter IJ so that the s-wavc interaction Aq is reduced by 
u = UN{0). The value of the parameter U is not equal 
to U, because the retardation and spin fluctuation effects 
should be taken into account. We consider [/ as a given 
parameter without estimating it microscopically. 

III. DEPENDENCE ON THE LAYER 
SPACING OF THE ANISOTROPIC PAIRING 
INTERACTION 

In this section, we calculate anisotropic components 
of the effective pairing interactions as functions of the 
layer spacing b. We concentrate on the case of intra-layer 
pairing ri{kz) = 1 for a while. 

The squared inverse of the screening length is 



9s = — p(m), 



(20) 



in Thomas-Fermi approximation, where p(/u) is total den- 
sity of states of electrons per unit volume at chemical po- 
tential II. In layer systems, the total density of states per 
unit volume p{fi) is written in terms of the total density 
of states per unit area (/i) in each layers as 



p(M) = pf(M)/6. 



(21) 



Here it is found that the screening becomes weaker when 
the layer spacing increases, because the volume density 
of electrons which contribute to screening decreases when 
the layer spacing increases. However, it should b(^ noted 
that the screening length within a layer changes by the 
change of the inter-layer spacing b, even when the lattice 
constant a in the layers is unchanged. Therefore the be- 
havior of the screening length examined is not derived by 
a simple scale transformation in terms of a and b as the 
length scales. 

We define a length scale bo as 



(22) 



from eqs. (20) and (21). In a simple case, the length scale 
bo is estimated as follows. Assuming non-interacting two 
dimensional electron gas in p':^^\ii), we obtain 



bo 



Trnan 



(23) 
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since p]^^\^) = mj-Kh^ and fcpa = ■\/27m, where n is the 
electron number per site. Here an denotes Bohr radius 
Oh = 47reo7i^/me^ = 0.5292A. As an example, if a '--^ 4A 
and n '-^ \ we have 60 9. 6 A as a crude estimation. 

Since the basic length scale an which is independent of 
the lattice constants a and h comes in eq. (23), changes 
not only of the ratio hja but also of the absolute values 
of a and h give rise to changes in the qualitative results. 




FIG. 1. The dimensionless coupling constants \m. as a 
function of the layer spacing 6. The solid and dashed lines 
show the results for p-wave (m = 1) and d-wave (m = 2), 
respectively. In the inset, the short dashed line shows the 
result for s-wave (m = 0). 
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to 1 in the limit of & = and decreases with h. It remains 
larger than the other anisotropic components, but if the 
additional short-range Coulomb energy U is sufficiently 
large so that Aq — u < Ai, p-wave pairing occurs instead 
of s-wave pairing. 

Figure ^ is the phase diagram at T = in h-U plane. 
It is found that p-wave superconductivity occurs in the 
region where the layer spacing 6 is larger and the short- 
range repulsion expressed by U is stronger. We will dis- 
cuss the reality of such parameter values in the layered 
compounds in the last section. 

On the other hand, for d-wave superconductivity to 
occur, some additional contribution to A2 or a negative 
contribution to Ai is needed, so that A2 becomes larger 
than Ai. We examine an enhancement of A2 due to an 
anisotropy of the density of states later, and briefly dis- 
cuss a contribution from the antiferromagnetic fluctua- 
tions in the last section. 

IV. INTER-LAYER PAIRING 

In this section we consider a situation in which the 
inter-layer coupling constant g' is of the same order as the 
intra-layer coupling constant g. The coupling constants 
would depend on the layer spacing 6, but here we regard 
them as independent parameters. The condition 5' ~ g 
would actually be satisfied when h is not too large. Then, 
we must consider the gap function of the form A(k) = 
A||(k||)77(fcz) with r\{\iz) = ^/2coskzb or ^/2smkzb. The 
expansion of A||(k||) by eq. (12) holds also in this case. 

Figure ^ shows the dimensionless coupling constants 
X'j^. A set of parameters, g' = O.Sg, IJ = OAg, and 
q'g = Qs are taken as an example. For b/bo ^ 0.6 and 
b/bo ^ 2.2, intra-layer pairing (of s-wave and p-wave in 
each region, respectively) is favored. On the other hand, 
for 0.6 ^ b/bo ^ 2.2, intcr-layer pairing with to = is 
favored. The gap function has a form such as 

J A(k) = Aosinfc^fe , . 

\ A(k) = Aq cos kzb. ^^^^ 

The formar is an order parameter of triplet pairing, while 
the latter is that of singlet pairing. These gap functions 
have horizontal line nodes at fc^ = 0, ±7r/6 and at fc^ = 
±7r/26, respectively, but they are isotropic in the layers. 



FIG. 2. The phase diagram at T = in b-O plane. SC 
stands for superconductivity. 

Figure |^ shows the result of X„i as a function of the 
layer spacing b. It is seen that both p-wave and d-wave 
components of the pairing interactions increase with the 
layer spacing b. In particular, it is found that the p-wave 
components increase rapidly in the region < 6 ^ &o- As 
the inset shows, the s-wave component Xo/gN{0) is equal 
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FIG. 3. The dimensionless coupling constants A^j of near- 
est neighbor layer pairing as a function of the layer spacing 
h. The thick solid and dashed lines show the results of the 
inter-layer pairing with m = and m = 1, respectively, while 
the thin dashed and short dashed lines show the results of 
intra-layer pairing, Ai and Aq — S, respectively. 



V. EFFECT OF ANISOTROPY IN THE 
ELECTRON DISPERSION 

In this section, we consider the square lattice systems, 
in which the electron dispersion depends on the direction 
of the momentum. We will show that the d-wave coupling 
constant A2 is enhanced for d,j.i_yi symmetry, but not 
for dxy symmetry, due to the anisotropy of the density 
of single-particle states. We define an angle-dependent 
density of states p(e, ip) as a density of single-particle 
states per unit energy and unit angle. 

In the square lattice system, the angle-dependent den- 
sity of states at the Fermi-energy p(0, (/s) can be approx- 
imated by 



p(0,V3) « po + P4 cos(4(/?), 



(25) 



where Lp is the angle between a momentum p and Px- 
axis. In addition, we regard fcp as being constant, for 
simplicity. 

Figure ^ shows a verification of this simplified model 
in the square lattice tight binding model with a nearest 
neighbor hopping energy t at /i = —t. Although the 
Fermi-surface is nearly isotropic, the density of states 
p(0, (/?) varies with the direction p. For example, when 
/i = —t, po ~ 0.142 and p4 ~ 0.040 are estimated. 

Regarding eq. (25) as an expansion of p(</3, 0), we 
could extend it into more general forms by adding terms 
P4„ cos(4rKp) with n > 2. Then the terms of p4„ mix 
of a small m with A^' of a large m! = |TO±4n|. However, 
since Vm decreases rapidly with m as seen by eq. (16), 
Am/s of such large vn! are small. Therefore the higher 
order terms in the expansion of p{jp^ 0) can be omitted in 
practice. 



In the gap equation, the anisotropic term propor- 
tional to p4 cos 4(y9 of the angle-dependent density of 
states mixes AmCosmi/j with A|„_4| cos[(to — 4)(^] and 
Am_|_4 cos[(r7i -f 4)(^], but it does not affect equations for 
AmSinTO(/3. Therefore we only consider equations for 
A,„. For general m, we can write the gap equation at 
T = as 



A^^^^og 
A„ + 

where we define A 



2e'>'wD 

P4 r Ur. 



-Am+4 + 



-^|m-4|} 



2po 'T'm+4 '^|m-4| 

(26) 

KnPo = KnA(O) is the dimen- 



(0) 



sionless coupling constant for the isotropic case. 




FIG. 4. The Fermi-surface of a tight binding model with 
a chemical potential fi = —t (solid line), and the averaged 
isotropic Fermi-surface with fcp « 1.97/a (short dashed line). 
The inset shows the angle-dependent density of states p(0, ip) 
at the Fermi-energy. Short dashed line in the inset shows the 
behavior of p(0, ip) approximated by eq. (25) with po ~ 0.142t 
and p4 = 0.040t. 

Since V3, V4, V5 • • • are much smaller than Vq and Vi, 
the terms proportional to p4 can be neglected in eq. (26) 
for m — and 1. Hence, Aq and Ai are not modified by 
P4. On the other hand, for m = 2, we cannot omit the 
term of A|„i_4| in eq. (26) since |m — 4| = 2. Neglecting 
the term of A,„+4 — Ag because Vg ^ V2 , we obtain 



Ao 



,(0) 



P4 

2po 



2e''tJD . 

log A2 



ttTc 



2e^wD . 
A2 log — ^ A2, 



(27) 



where we define an effective coupling constant A2 = 
A^°^(l + P4/2po), which gives by eq. (17). 

Therefore, it is found that c?a;2„j,2-wave pairing is fa- 
vored more than dxy-'wave pairing by the enhancement 
factor (l + p4/2po), since p4 changes the equation for A^ 
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but does not for as mentioned above. The enhance- 
ment factor l-|-p4/2po is estimated to be 1.14 for fj, = —t, 
and 1.22 for ^ = —0.5t. On the competition with p-wave 
pairing, those values are not large enough to change the 
sign of A2 — Ai. Therefore, another non-phonon contri- 
bution seems to be needed for d-wave pairing to occur. 

Figure |^ are phase diagram in the absence of p-wave 
superconductivity. It is found that d-wave superconduc- 
tivity is favored when the layer spacing b is larger and 
the short-range Coulomb repulsion U is stronger. If there 
are additional contributions to d-wave superconductivity 
mentioned above, the phase boundary is shifted down- 
ward. 



J-wave SC 



s-wave SC 
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b/bo 



The phase diagram of s-wave superconductivity 
j^2-wave superconductivity at T = in b-U plane 



FIG. 5 

and d^2_ 

in the absence of p-wave superconductivity. SC stands for 
superconductivity. 



VI. DISCUSSION AND SUMMARY 

We have examined pairing interactions mediated by 
phonons in the layer systems. The screening of Coulomb 
interactions becomes weaker when the layer spacing b 
increases. Then anisotropic components of the pairing 
interactions increase with the layer spacing b since the 
momentum dependence of the interactions changes. In 
particular, p-wave superconductivity occurs for large b 
and strong short-range Coulomb repulsion [/, even in the 
absence of any additional non-phonon interactions. 

It was found that the p-wave coupling constant Ai in- 
creases rapidly with the layer spacing b in the region 
b ^ bo, where bo is a length scale defined by eq. (22). 
For the rapid increase of Ai, the condition Aq — u < Ai is 
realized more easily in layer systems than in usual three 
dimensional systems, where Aq denotes the s-wave cou- 
pling constant and ~u is a negative contribution to s- 
wave pairing due to the short-range Coulomb repulsion 
discussed near eq. (19). Hence triplet pairing supercon- 
ductivity is favored in layered compounds. 



We have also examined inter-layer pairing. In some re- 
gion of the parameter space, for example 0.6 ^ b/bo ^ 2.2 
for the parameters indicated in Fig. ^, the gap function 
has horizontal line nodes parallel to the layers. In this 
case, the solutions of singlet pairing and triplet pairing 
of eq. (24) degenerate. If some effect due to spin fluctu- 
ations, ferromagnetic correlations, or spin-orbit coupling 
removes this degeneracy, inter-layer triplet pairing may 
occur. 

In Sr2Ru04 compounds, existence of the line nodes 
was supported by some experiments such as tempera- 
ture dependences of specific heat and NMR relaxation 
rate |27|. However, direction of the line nodes does 
not seem clear at the present. Line nodes vertical to 
the layers were indicated by ultrasound attenuation , 
whereas almost isotropic state was indicated by thermal 
conductivity [^. The isotropic state can be consistent 
with the specific heat and NMR experiments, if the hor- 
izontal line nodes are assumed 

In the present theory, inter-layer triplet pairing re- 
produces the horizontal line nodes, while the intra-layer 
triplet pairing is a candidate for the vertical line nodes. 
For the latter pairing, we need some additional mecha- 
nism for the vertical line nodes to occur, because isotropic 
states such as px + ipy have the lowest free energy in the 
present isotropic system. Consistent explanation of the 
experimental results within the present theory remains 
for a future study. 

In order to discuss the reality of the phonon-mediated 
anisotropic superconductivity, we crudely estimate the 
parameters for the Sr2Ru04 compound and quasi-one- 
dimensional organic superconductors from the observed 
transition temperature Tc ~ 1.5K. We assume triplet 
pairing here, although for the organics it might be rather 
controversial. The results of the parameter values do not 
strongly depend on the direction of the line nodes. 

Roughly speaking, b > bo is satisfied in both kinds of 
compounds. Thus we have Ai ~ 0.15 x gN{0) for intra- 
layer triplet pairing from Fig. |l|, while Aq ^ 0.22 x gN{0) 
for inter-layer pairing from Fig. |^. On the other hand, 
if we assume wd ~ lOOOK and Tc ~ 1.5K, we have 
Ai ~ 0.151 (or Aq ~ 0.151). Therefore, we obtain 
gN{0) ~ 1.0 and 0.69, respectively. For such choices of 
parameter values, it is likely that s-wave pairing is sup- 
pressed. For example, for intra-layer pairing, since we 
obtained gN{0) ~ 1 above, the s-wave coupling constant 
is estimated as Aq ~ 0.58 x gN{0) ^ 0.58 at 6 = 5o by 
eq. (16). Thus the difference of Aq and Ai is about 0.42, 
which corresponds to Vo - Vi - 0.42/7V(0) - 0.421^, 
where W is the band width. Therefore s-wave pairing 
is suppressed when the magnitude of the negative con- 
tribution U due to the short-range Coulomb repulsion is 
larger than 0A2W ~ W/2. Although this estimation is 
crude, the value ~ W/2 seems realistic as the order of 
the magnitude. 

On the other hand, for d-wave superconductivity to 
occur in the present model, it is needed that p-wave pair- 
ing is suppressed for some extra reason as well as s-wave 
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pairing, or that there are some additional contributions 
to d-wave pairing. For this problem, we examined effect 
of the anisotropy of the electron dispersion. It was found 
that the d-wave coupling constant A2 is enhanced by the 
anisotropy for dj.2_y2 symmetry, while not for dxy sym- 
metry and Px, Py symmetries. However, the enhancement 
does not seem to be large enough to realize the d-wave 
superconductivity. 

This might suggest an existence of a non-phonon 
contribution to the c?-wave pairing interaction in the 
cuprates. For example, many authors discussed that 
antiferromagnetic spin fluctuations may contribute to d- 
wave pairing interactions pl| , p^ ,p|,p| j3^j20| , ^ . In partic- 
ular, such interactions enhance d-wave pairing efflciently 
in the presence of the van Hove singularities in the 
square lattice systems. 

However, even if we assume that a non-phonon con- 
tribution is indispensable for high-Tc, the present theory 
suggests that there is a large phonon contribution to the 
d-wave pairing interactions especially in layer systems for 
the weak screening. This result is consistent with the 
observed large shifts of as absolute values due to the 
isotope effect ||To|-[l6t. 

It was also found that the coupling constant A2 
increases with the layer spacing h. This behav- 
ior might be a reason why the transition tempera- 
ture of Bi2Sr2CaCu208+5 is much higher than that of 
La2-2;Sr3;Cu04. Since is a sensitive function of A,„, 
such a slight enhancement of A2 may increase consid- 
erably. 

In conclusion, phonon-mediated pairing interactions 
include anisotropic components of various symmetries. 
In particular, they are large when the screening effect is 
weak due to the layered structure of the system. As a re- 
sult, the interactions induce various types of anisotropic 
superconductivity depending on values of the energy pa- 
rameters and lattice constants. In particular, triplet su- 
perconductivity is favored for large layer spacing and 
strong short-range repulsion. For anisotropic supercon- 
ductivity, Tc increases with the layer spacing 6. The 
phonon-mediated pairing interactions may play some es- 
sential roles in the anisotropic superconductivity of lay- 
ered compounds, such as Sr2Ru04, the organic supercon- 
ductors, and the high-Tc cuprates. 
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